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ABSTRACT  OF  FINAL  REPORT 


During  the  course  of  this  research  contract,  we  have  concentrated  on  the  device  and  system  aspect 
of  distributed  antenna  systems  aimed  towards  sipport  wireless  communications  in  an  in-building 
environment  The  infrastructure  is  fiber-based,  but  with  the  key  criterion  being  low-cost.  As  such, 
we  address  the  issue  both  from  the  device  end  and  from  the  system  end.  From  the  device  aid,  we 
investigated  VCSELS  as  optical  sources  for  the  fiber-optic  network,  since  VCSELS  are  the  lowest 
cost  of  ANY  active  transmitter  component  available  to  date,  and  projected  to  continue  into  the 
fiiture.  In  particular,  we  investigated  VCSELS  novel  structures  in  which  a  voltage  biased  quantum 
well  absorber  is  embedded  in  the  mirror  stack.  The  added  RF  functionalities  lead  to  several 
appUcations  for  these  devices  were  demonstrated. 

From  the  systems  end,  we  have  examined  network  architectures  which  lead  to  the  lowest  cost 
solution.  The  architecture  we  developed  is  called  a  “saturated”  antenna  architecture.  This  approach 
offers  the  possibility  of  negating  the  need  for  site-specific  planning,  modeling  and  survey  for 
dqiloyment  of  in-building  wireless  S5^ems,  and  is  possible  only  with  RF  photonics  hardware  of 
sufficifflitly  low  cost,  such  as  those  offered  by  VCSEL  transmitter  systems.  We  investigate  the 
characteristic  radio  distribution  statistics  under  this  scenario,  and  found  that  with  a  "sufficiently 
large"  number  of  antennae  elements  per  unit  area,  the  radio  statistics  approaches  certain  limiting 
distributions  which  are  largely  ind^aidoit  of  the  specific  flcxjr-plan.  These  limiting  distributions 
(sn  be  accuratdy  predicted,  without  invokiiig  adjustable  parameters,  using  only  a  few  key 
"global"  building  parameters  which  can  be  easily  determined  for  a  given  buildiDg.  A  system  of  this 
type  can  thus  be  rapidly  and  inexpensively  dqjloyed. 

Section  1  of  this  rqjort  describes  in  detail  the  design  and  RF  characteristics  of  cjuantum-well 
absorber-embedded  VCSELS,  and  the  RF  photonic  appheations  made  possible  by  this  novel 
device.  Sechcjn  2  describes  the  basic  “saturated”  antama  architecture  approach,  and  extensive 
simulation  results  justifying  the  high  performance  of  this  approach.  Section  3  approaches  the 
“saturated”  antama  architecture  analytically,  and  arrives  at  a  fundamaital  imderstanding  of  the 
mechanism  on  which  the  saturated  antama  network  operates,  and  factors  which  leads  to  the  results 
we  observed  in  the  simulaticm.  This  approach  leads  to  a  sohd  unfestanding  of  the  performance  of 
such  a  network  and  its  dependence  cm  various  parameters  of  the  system. 
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1.  RF  photonics  by  VCSELS  with  embedded  quantum-well  absorber 
Section  Abstract 

A  VCSEL  with  an  embedded  quantum-well  absorber  exhibits  very  interesting  and  useful  RF 
properties.  The  underpinning  of  these  properties  are  regimes  of  n^ative  differential  resistance  in 
the  absorber  under  varying  laser  bias  conditions.  By  choosing  the  absorber  bias  conditions  such 
that  the  operating  load  line  intersects  the  absorber  IV  trace  three  times,  the  laser  can  be  made  to 
exhibit  optical  bistabihty.  The  magnitude  of  the  bistability  and  consequent  hysteresis  loop  can  be 
controlled  by  adjusting  the  absorber  bias.  If  the  bias  conditions  are  adjusted  so  that  the  load  line  is 
tangmtial  to  the  portion  of  the  absorber  I-V  exhibiting  native  differential  resistance,  then  the 
device  can  be  made  to  self-pulsate.  Self-pulsations  at  frequencies  as  hi^  as  2  GHz  with  700  kHz 
RF  linewidth  FWHM  (full-width  at  half-maximum)  have  been  obtained  using  these  devices. 
Furthermore,  the  self-pulsation  frequency  can  be  tuned  over  700  MHz  by  adjusting  the  bias 
conditions,  representing  a  substantial  advance  over  existing  self-pulsating  VCSELs. 

We  have  also  demonstrated  a  novel  modulation  sdieme  using  these  devices,  in  vduch  the  drive 
signal  is  applied  to  the  bias  voltage  across  the  absorber.  Theoretical  analysis  of  the  chirping 
mechanism  leads  us  to  expect  that  this  technique  will  minimize  the  chirp  at  hi^  modulation 
frequencies,  while  still  providing  substantial  modulation  depth  and  hi^  speeds.  We  achieved  a  -3 
dB  bandwidth  of  9  GHz;  the  bandwidth  increases  with  the  laser  bias  current  at  a  fate  of  7 
GHz/mA^°. 

1.1  Introduction 

We  have  constructed  a  novel  vertical-cavity  surface-emitting  laser  (VCSEL)  with  a  quantum-well 
absorber  integrated  into  the  mirror  stack.  The  absorpticm  can  be  varied  by  applying  a  voltage  bias 
across  the  absorber  junction,  which  effectively  changes  the  mirror  reflectivity  and  hence  varies  the 
laser  output  power.  We  have  demonstrated  several  apphcatioos  for  this  new  device,  including 
optical  bistability,  self-pulsation,  and  a  novel  intracavity  absorber  modulation  scheme. 

Self-pulsating  lasers  have  found  extensive  commercial  use  in  applications  such  as  optical  disk 
readout  and  multimode  fiber  communications  networks  because  the  extremely  short  coherence 
length  of  self-pulsating  lasers  renders  them  relatively  insensitive  to  optical  feedback.  Althou^ 


2 


self-pulsation  has  been  demonstrated  previously  in  edge-emitting  lasers  which  use  int^ated  active 
saturable  absorbers  to  obtain  controllable  self-pulsations  up  to  GHz  frequencies,  recent  work  has 
focused  on  obtaining  self-pulsating  VCSELs.[l-3]  The  surface-nonnal  geometry  of  VCSELs  leads 
to  advantages  such  as  circular  beam  output  and  the  potential  for  two-dimensional  arrays.  Initial 
demonstrations  of  self-pulsating  VCSELs  used  defects  caused  by  ion  implantation  as  passive 
saturable  absorbers;  these  efforts  produced  at  most  MHz-range  pulsations  with  RF  Unewidths  of 
rou^y  100  MHz.  [4]  Using  our  new  VCSEL,  which  has  an  active  quantum-well  absorber 
int^ated  into  the  mirror  stack,  we  obtain  reproducible  self-pulsations  with  repetition  rates  of 
greater  than  2  GHz  and  RF  linewidtlis  of  700  kHz.  [5,6]  Furthermore,  we  demonstrate  tuning  of 
the  self-pulsation  frequency  over  700  MHz.  In  this  section,  we  present  the  experimental  results 
and  discuss  in  detail  the  mechanism  which  yields  optical  bistability  and  self-pulsation  in  these 
devices. 

We  also  use  this  structure  to  demonstrate  a  novel  modulation  technique  in  which  the  modulation 
signal  is  applied  directly  to  the  intracavity  quantum-well  absorber.  Ideally,  a  VCSEL  modulation 
scheme  should  provide  a  hi^  modulation  depth  and  low  frequency  chirp,  while  still  being  hi^- 
speed.  Although  direct  current  modulation  has  been  demonstrated  at  hi^  speeds  and  for  large 
modulation  depths,  changing  the  carrier  density  results  in  frequency  chirp,  limiting  both  digital  and 
large-signal  AM  applications.  [7-9]  Frequency  chirp  can  be  avoided  by  using  an  external 
modulator;  however,  the  extinction  ratio  obtained  using  this  technique  is  low  and  insertion  loss  is 
problematic.  As  an  alternative,  we  demonstrate  a  scheme  in  which  the  modulation  is  apphed  to  a 
quantum-well  absorber  integrated  the  VCSEL  mirror  stack.  [10-1 1]  Modulating  the  absorber 
causes  the  mirror  reflectivity  to  vary,  resulting  in  modulation  of  the  laser  output  power.  We  expect 
the  modulation  bandwidth  using  this  technique  to  be  comparable  to  that  of  direct  current 
modulation.  The  device  can  also  be  designed  for  hi^  modulation  depth  by  positioning  the 
modulator  quantum-well  absorption  edge  to  be  shghtly  shorter  than  the  Fabry-Perot  wavelength,  so 
that  as  the  absorption  edge  red-shifts  with  increasing  apphed  reverse  bias  it  is  swept  across  the 
Fabry-Perot  wavelength.  The  placement  of  the  absorber  at  a  peak  of  the  internal  optical  field 
intensity  distribution  results  in  high  modulation  efficiency.  Furthermore,  we  analyze  the 
modulation  chirp  properties  and  show  that  our  technique  should  lead  to  significantly  less  frequency 
chirp  than  does  direct  modulation  at  hi^  frequencies.  Using  intracavity  absorber  modulation,  we 
demonstrate  a  3-dB  bandwidth  of  9  GHz  and  a  modulation  efficiency  of  7  GHz/VmA.  These 
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results  are  liinited  only  by  the  relaxation  osotiation  frequency  and  could  be  increased  substantially 
by  inqjlemoiting  the  design  in^rovements  suggested  below. 

1.2  Device  Structure 

The  device,  which  is  shown  in  Figure  la,  is  a  VCSEL  with  an  n-p-n  configuration.  The  bottom  p-n 
junction  is  used  to  forward  bias  a  two  or  Ihree  quantum-well  gain  region.  The  top  n-p  junction 
allows  reverse  biasing  of  a  sin^e  quantum-well  absorber  positioned  in  the  top  mirror  stack  as 
shown.  Current  and  modal  confinement  are  achieved  using  two  oxide  layers.  This  device  is  similar 
to  that  of  an  int^rated  VCSEL  detector  in  Reference  12,  where  details  of  design  and  fabrication 
can  be  found.  [12] 


The  electronic  driving  circuit  used  with  the  device  varies  depending  on  the  desired  application.  To 
obtain  optical  bistability  or  self-pulsation  of  the  device,  the  bias  circuit  shown  in  Figure  lb  is  used. 
The  absorber  is  biased  using  a  variable  resistor  R  in  series  with  a  voltage  source  Vo.  A  50Q 
termination  path  is  provided  for  any  AC  current  generated  in  the  absorber,  and  the  gain  section  is 
biased  with  a  DC  current  source.  The  driving  circuit  used  to  perform  hi^  speed  absorber 
modulation  is  shown  in  Figure  Ic.  In  this  case,  the  absorber  is  reverse  biased  with  a  DC  voltage 
source,  and  the  AC  modulation  signal  is  siqierinposed  using  a  bias-T.  The  p-contact  of  the  device 
serves  as  the  common  ground  for  both  the  absorber  and  the  active  r^ons.  A  DC  current  source  is 
used  to  punp  the  active  region,  wbile  a  50-fi  termination  path  is  provided  for  any  AC  signal 
generated  in  the  gain  region. 


n-DBRs 
QW  absorber 
oxide  layer 

2  or  3-QW 
active  region 

oxide  layer 


Figure  la.  Device  schematic 
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Figure  lb.  Electronic  bias  circuit  for  optical 
bistability  and  self-pulsation. 


50  Q 


^aAC  ^a,DC 


Figure  Ic.  Electronic  bias  circuit  for  absorber 
modulation. 


1.3  Negative  Differential  Resistance 

The  longitudinal  modes  of  this  device  are  very  widely  spaced  due  to  the  diaracteristic  short  cavity 
length  of  a  VCSEL.  Thus,  the  lasing  wavelength  of  the  device  at  a  fixed  bias  current  is  determined 
by  the  Fabry-Perot  wavelength,  providing  a  design  flexibihty  which  is  not  available  in  edge- 
emitting  devices.  This  flexibility,  which  is  discussed  in  Reference  6,  allows  us  to  obtain  a  large 
negative  differential  resistance  in  the  absorber  using  small  absorber  bias  voltages  and  potentially  at 
hi^  laser  powers[6]. 

To  illustrate  this  effect,  we  engineered  two  VCSELs,  with  differing  placements  of  the  lasing 
wavelength  relative  to  the  absorption  wavelength,  defined  by  the  peak  of  the  absorption  spectrum 
when  no  bias  voltage  is  applied  to  the  absorber.  In  Device  A,  the  absorption  wavelength,  which  is 
designed  to  be  851nm,  is  23nm  longer  than  the  lasing  wavelength  of  828nm.  The  gain  r^on  of 
this  device  is  composed  of  three  70A  Ino  2Gao.gAs  quantum-weUs  and  is  designed  to  be  centered  at 
843nm.  The  single  quantum-well  absorber  is  90A  thick,  with  the  current-voltage  (k-Va) 
characteristic  shown  in  Figure  2a.  In  Device  A,  the  la-Va  response  shifts  downwards  with 
increasing  laser  power,  so  the  absorber  behaves  like  a  typical  photodiode. 
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In  contrast.  Device  B  is  designed  so  that  the  absorption  wavelength  (approximately  940ran)  is 
26nm  shorter  than  the  lasing  wavelength  (966nm).  The  gain  r^on  centered  at  940nm  is 
comprised  of  two  80A  quantum-wells,  and  the  absorber  quantum  well  is  also  80A  thick.  The 
absorber  T-Va  characteristic  for  Device  B  exhibits  a  regime  of  n^ative  difiPerential  resistance,  as 
shown  in  Figure  2b.  Figure  3  illustrates  the  origin  of  this  regime.  For  a  constant  DC  gain  bias 
current,  the  lasing  wavelength  is  fixed.  When  a  voltage  bias  is  applied  to  the  absorber,  the 
absorption  spectrum  reddiifts  and  broadens  due  to  the  quantum  confined  Stark  effect,  effectively 
sweeping  the  absorption  peak  across  the  lasing  wavelength  as  shown  in  Figure  3a.  This  causes  a 
r^on  of  n^ative  differential  resistance  in  the  absorber  k-Va  response;  see  Figure  3b.  Typically, 
an  absorber  bias  voltage  (Va)  of  about  1.5V  is  required  to  obtain  the  negative  differential 
resistance.  This  bias  voltage  is  appHed  across  the  5X/4  (300nm)  intrinsic  r^on,  so  the  apphed 
field  is  about  5x10"^  V/cm. 


Furthermore,  the  lasing  wavelength  redshifts  with  applied  laser  bias,  altering  its  position  relative  to 
the  absorption  wavelength.  This  observation  suggests  the  potential  to  obtain  high  power  operation 
of  our  device  by  positioning  the  Fabry-Perot  wavelength  sU^tly  closer  to  the  absorption  peak. 


Figure  2a.  Absorber  current  vs.  voltage  (k- 
Va)  traces  for  Device  A,  with  varying  laser 
bias  current  k-  device,  the 

absorption  wavelength  is  23nm  longer  than 
the  lasing  wavelength. 


absoiber  voltage  \  (V) 

Figure  2b.  Absorber  current  vs.  voltage  (k- 
Va)  traces  for  Device  B,  with  varying  laser 
bias  current  k-  For  this  device,  the 
absorption  wavelength  is  26nin  shorter  than 
the  lasing  wavelength.  The  k-Va  traces 
exhibit  a  regime  of  sharply  negative 
differential  resistance  (NDR). 
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Figure  3a.  Origin  of  the  sharp  region  of 
negative  differential  resistance  in  the  absorber  I- 
V  trace.  The  absorption  spectrum  redshifts  with 
applied  voltage  bias,  sweeping  the  absorption 
peak  across  the  lasing  wavelength. 


Figure  3b.  The  absorber  IV  trace  resulting 
from  the  sweeping  of  the  absorption  peak  in  Fig 
3a  displays  negative  differential  resistance. 


1.4  Optical  Bistability 

We  exploit  the  negative  differential  resistance  exhibited  by  the  absorber  la-Va  traces  of  Device  B  in 
order  to  obtain  either  optical  bistability  or  self-pulsation.  For  these  apphcations.  Device  B  is 
electronically  biased  as  shown  in  Figure  lb.  Our  work  on  optical  bistabihty  is  presented  in 
Reference  5. [5] 

The  device  can  be  operated  in  a  bistable  manner  by  choosing  the  absorber  bias  (R  and  Vo)  such 
that  the  operating  load  line  intersects  a  single  absorber  k-Va  trace  three  times,  as  shown  in  Figure 
4a. [3,5]  Power-current  (L-Ig)  traces  of  a  typical  Device  B,  with  the  absorber  biased  so  that  the 
load  line  intersects  the  absorber  k-Va  trace  three  times,  are  shown  in  Figure  4b.  A  noticeable 
discontinuity  in  dL/dIg  and  a  hysteresis  loop  occur  at  threshold,  indicating  bistable  operation. 
Figure  4b  shows  that  the  hysteresis  loop  can  be  closed  by  decreasing  the  absorber  bias  resistance 
R. 
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Figure  4a.  Bistable  operation;  the  absorber  is 
biased  so  that  the  operating  load  line  intersects 
the  absorber  la-Vg  trace  three  times. 


Figure  4b.  Bistable  operation:  the  L-lg  trace 
exhibits  a  discontinuity  in  dL/dIg  and  a 
hysteresis  loop.  Absorber  is  biased  with 
Vo=40V  and  R=500kn,  350kQ,  or  300kf2. 


1.5  Self-pulsation 


Self-pulsation  of  Device  B  is  achieved  by  adjusting  the  absorber  bias  conditions  so  that  the 
operating  load  line  is  tangential  to  the  region  of  negative  differential  resistance  in  the  absorber  la-V, 
trace;  see  Figure  5a.  Because  Device  B  is  designed  as  discussed  in  Section  3  to  obtain  r^ons  of 
sharp  native  differential  resistance  in  the  absorber  h-V,  characteristic,  relatively  low  absorber 
biases  (R=10-100kf2,  Vo=5-15V)  are  required  to  obtain  self-pulsation.  A  complete  discussion  of 
our  work  on  obtaining  self-pulsation  is  given  in  References  5  and  6.  [5,6] 


Figure  5a.  Self-pulsation  regime:  the  absorber  is 
biased  so  that  the  operating  load  line  is 
tangential  to  the  absorber  la-Va  trace. 


gain  current  Ig  (mA) 


Figure  5b.  Self-pulsation  regime:  the  L-lg  trace. 
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Typical  power  versus  current  (L-Ig)  traces  are  shown  in  Figure  5b  both  for  the  case  in  which  the 
absorber  bias  is  turned  off  (R=0,  Vo=0),  so  that  the  device  operates  as  a  standard  VCSEL,  and  also 
for  the  case  in  which  the  absorber  is  biased  appropriately  for  self-pulsation.  As  expected,  when  the 
absorber  bias  is  turned  on,  the  lasing  threshold  of  the  device  increases.  For  the  particular  absorber 
bias  (R,  Vo)  used,  the  L-Ig  curve  exhibits  a  sharp  jump  and  a  narrow  hysteresis  loop  near  threshold. 
We  observe  self-pulsation  of  the  device  when  the  laser  bias  current  is  about  1 .5  times  the  threshold 
current,  as  indicated  by  the  boxed  region  in  Figure  5b.  The  jump  in  the  L-lg  curve  indicates  that 
just  above  threshold  the  device  is  operating  in  a  bistable  manner  because  the  operating  load  line 
thrice  intersects  the  absorber  h-Va  trace  when  the  gain  bias  current  Ig=5mA.  When  the  gain 
current  lg=6-7mA,  which  is  the  necessary  condition  for  self-pulsation  with  the  given  absorber  bias 
conditions  (R=75kn,  Vo=7V),  the  load  line  is  tangaitial  to  the  absorber  la-V,  response.  In 
practice,  self-pulsation  of  the  device  can  be  obtained  at  a  range  of  bias  currents  (Ig)  and  output 
powers  by  choosing  R  and  Vo  so  that  the  load  line  is  tangential  to  the  absorber  Ij-V,  trace  at  the 
desired  bias  current.  This  unique  feature  arises  due  to  the  abUity  to  control  the  laser  and  absorber 
wavelengths  independently,  rqpresenting  a  major  design  strength  of  our  device. 

We  measure  the  output  of  the  device  using  an  RF  spectrum  analyzer.  When  the  device  is  biased  so 
that  it  does  not  self-pulsate,  the  RF  spectrum  show/s  only  a  broad  peak  due  to  relative  intensity 
noise;  see  Figure  6a.  When  the  absorber  bias  conditions  are  adjusted  such  that  the  operating  load 
line  is  tangential  to  the  region  of  native  differential  resistance  in  the  absorber  IV,  the  device  self- 
pulsates  at  a  GHz-range  repetition  rate,  resulting  in  the  sharp  peak  in  the  RF  spectrum  shown  in 
Figure  6b.  The  self-pulsation  peak  is  rou^y  40dB  above  background,  and  has  an  extremely 
narrow  RF  linewidth  of  700  kHz,  as  illustrated  by  the  enlarged  RF  trace  shown  in  Figure  6c. 

The  observed  self-pulsations  occur  at  a  repetition  rate  equal  to  the  relaxation  oscillation  frequency 
of  the  device.  Hence,  the  self-pulsation  frequaicy  can  be  tuned  by  varying  the  DC  gain  bias 
current.  As  the  gain  bias  current  is  varied,  the  absorber  bias  (R  and  Vo)  must  be  adjusted  so  that 
the  operating  load  line  remains  tangaitial  to  the  absorber  IV  in  order  to  maintain  self-pulsation. 
Figure  7  shovws  that  the  pulsation  frequency  obeys  a  square  root  dependence  on  the  gain  bias  above 
threshold,  as  expected.  We  demonstrate  tuning  of  the  self-pulsation  frequency  over  700  MHz 
using  this  technique. 
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frequency  (GHz) 

Figure  6a.  Self-Pulsation  OFF;  RF  Spectrum 
Analyzer  trace  of  the  device  output  when  the 
absorber  is  not  biased  for  self  pulsation 
(R=100kQ,  Vo=10V,  Ig=5.3mA). 


).4  1.4  2.4 

frequency  (GHz) 

Figure  6b.  Self-Pulsation  ON:  RF  Spectrum 
Analyzer  trace  of  the  device  output  when  the 
absorber  is  biased  for  self  pulsation  (R=75kO, 
Vo=7V,  Ig=6mA). 


frequency  (GHz) 

Figure  6c.  Self-Pulsation  ON;  RF  Spectrum  enlarged  to  show  RF  linewidth  of 
self-pulsation  peak. 


Figure  7.  Tuning  of  the  self-pulsation  frequency  by  adjusting  the  gain  bias 
current  Ig. 
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1.6  absorber  modulation 


1.6.1  Experimental  Results 

We  demonstrate  hi^-speed  modulation  of  the  intracavity  absorber  using  a  device  of  Type  A  (see 
Figure  2a  above).  This  work  is  described  in  References  10  and  11.  [10, 11]  The  828nm  lasing 
wavelength  of  Device  A  is  23nm  shorter  than  the  absorption  wavelength,  which  is  designed  to  be 
851nm.  The  electronic  driving  circuit  used  for  these  measurements  is  shown  in  Figure  Ic. 

The  L-Ig  and  Vg-I^  charactenstics  of  a  typical  Device  A  are  shown  in  Fig.  8.  The  threshold  current 
is  1 .7  mA,  and  the  threshold  voltage  is  5.6  V.  The  output  power  typically  peaks  at  0.7  mW. 


Figure  8.  L-Ig  and  Vg-Ig  for  a  typical  Device  A 
with  the  absorber  voltage  Va=0. 


Figure  9a  shows  the  response  of  the  device  to  absorber  modulation  under  several  DC  bias 
conditions  (Ig)  for  the  gain  r^on.  The  experimental  results  indicate  a  -3-dB  bandwidth  of  9  GHz 
at  Ig  =  8  mA.  The  bandwidth  as  well  as  the  relaxation  oscillation  frequaicy  increase  with  laser 
bias  current  Ig  as  expected  in  conventional  modulation;  see  Figure  9b.  We  achieved  7  GHz/VmA  in 
the  linear  r^me,  rolling  off  to  1 .7  GHz/VmA. 
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Figure  9a.  Absorber  modulation  response  for 
fixed  DC  absorber  bias  (Vajxf="l  V)  and  varying 
laser  bias  current 
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Figure  9b.  L-Ig  and  V^-Ig  for  a  typical  Device  A 
with  the  absorber  voltage  V^=0. 


we  also  measured  the  absorber  modulation  response  for  fixed  laser  bias  currents  and  at  different 
DC  absorber  biases  The  results  iir  Figure  10  indrcate  that  the  bandwidth  does  not  change 
appreciably  with  varying  absorber  DC  bias. 


frequency  (GHz) 


Figure  10.  Absorber  modulation  response  for 
fixed  gain  bias  current  Ig  and  varying  DC 
absorber  bias  (Vaxc  =  OV,  -IV,  -3V,  -4V). 


1.6.2  Discussion  of  Frequency  CMrp 

The  ftejuency  chirp  due  to  intracavity  absorber  modulation  is  compared  in  the  following  discussion 
and  in  Reference  11  to  the  chirp  introduced  by  conventional  diiect  current  modulation.[lll 
Modulating  the  absorber  effectively  modulates  the  mirror  loss  since  the  absorber  is  part  of  the 
VCSEL  mirror  stack.  This  in  turn  leads  to  modulation  of  the  photon  lifetime. 
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We  compare  the  ftecprency  chirp  due  to  intracavity  absorber  modulation  (dip,  vdiere  Xp  is  the  photon 
lifetime)  to  the  chirp  obtained  using  direct  current  modulation  (5J,  where  J  is  the  injected  current 
density)  by  performing  a  small-signal  analysis  of  the  rate  equations.  Define  Sn^  to  be  the  change  in 
the  carrier  density  of  the  gain  region  when  the  absorber  is  modulated  (i.e.  5J  =  0).  Likewise,  define 
5nj  to  be  the  change  in  the  carrier  density  of  the  gain  region  when  the  pump  current  is  modulated 
(i.e.  6xp  =  0).  The  ratio  of  the  two,  |  Sn^  /  5nj  |,  provides  a  comparison  of  the  chirp  produced  by  the 
two  techniques.  Note  that  this  ratio  depends  on  the  modulation  frequency  o; 


f  dg  ^ 

Sv  r^  +  2r^,^Bn 

\  on  ) 

f .  1  „  dg^ 

IG)  -|-  -f  gVa 

V  lAe 

**  Tp  *  asj 

(1) 

where  Vg  =  group  velocity,  S  =  photon  density,  g  =  gain,  F  =  confinement  factor,  Psp  =  spontaneous 
emission  factor,  x^e  =  (A+2Bn+3Cn^)'',  A  =  nonradiative  recombination  coefficiait,  B  =  radiative 
recombination  coefficient,  and  C  =  Auger  recombination  coefficient. 


The  magnitude  of  the  two-pole  response  given  by  (1)  monotonically  decreases  with  or  due  to  the 
heavy  damping  introduced  by  5g/5S  and  5g/dn.  For  frequencies  greater  than  fc  (the  frequency  at 
which  I  5n./5nj  |  =  1),  absorber  modulation  produces  less  chirp  than  direct  modulation.  From 
equation  (1): 
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where  cor  is  the  resonance  frequency.  Since  5g/5S  is  negative,  fc  is  less  than  the  resonance 
frequency.  Figure  1 1  shows  the  calculated  effect  of  absorber  and  gain  modulation  on  the  carrier 
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density  in  the  gain  r^on  for  the  same  photon  modulation  depth  and  standard  VCSEL  parameters. 
From  Fig.  13,  the  ratio  of  the  -3dB-frequency  f.3dB  to  fc  is  1.6,  so  the  critical  frequency  occurs  well 
within  the  useful  bandwidth  of  the  device.  Hence,  absorber  modulation  at  frequencies  between  fc 
and  f.3dB  results  in  less  chirp  than  does  direct  gain  modulation.  The  carrier  density  within  the 
absorber  layer  is  also  modulated.  However,  the  absorber  is  far  away  from  the  active  r^on,  so  any 
frequeicy  chirping  effects  should  be  quite  small. 


Figure  11.  Predicted  effect  of  absorber  and  direct  modulation  on  carrier  density  in 
the  gain  region.  When  |5ny|<|5nj|,  absorber  modulation  results  in  less  chirp  than  does 
direct  gain  modulation. 


1.6.3  Design  Criteria 

The  relative  placement  of  the  gain  peak  wavelength  (Again),  the  absorption  peak  wavelength  (Aabs), 
and  the  Fabry-Perot  transmission  peak  wavelength  (Aep)  is  an  important  design  criterion  for  this 
device.  Ideally,  we  should  have  Again  App  in  order  to  achieve  maximum  overlap  of  the  gain  peak 
with  the  Fabry-Perot  as  the  gain  peak  red-shifts  with  increasing  pump  current.  Furthermore,  we 
desire  that  Aabs  <  App  with  these  two  wavelengths  fairly  close  to  each  other,  so  that  a  small  diange  in 
the  apphed  voltage  across  the  absorber  induces  a  large  change  in  the  absorption  due  to  the  quantum- 
confined  Stark  effect.  Since  the  wavelength  of  the  pumped  gain  region  will  shift  ftuther  than  will  the 
wavelength  of  the  voltage-biased  absorber,  the  ideal  relative  positioning  of  the  three  wavelengths 
should  be  Agai„  <  Aabs  <  App. 
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Next  consider  vdiat  happens  v\hm  the  ideal  conditions  are  not  met.  Placing  the  Fabry-Perot 
wavelaigth  too  much  longer  than  the  peak  wavelength  of  the  gain  or  absorber  leads  to  poor 
modulation  efficiency,  as  shown  in  Reference  14,  since  greater  voltages  across  the  absorber  are 
required  to  swe^  the  absorption  edge  across  the  Fabry-Perot  wavelength.  [14]  Furthermore,  the 
smaller  overlap  of  the  gain  peak  with  the  Fabry-Perot  results  in  a  narrow  bandwidth,  also  evident  in 
Reference  14. 

In  the  device  used  for  the  experimental  measuremaits  discussed  in  Section  6.1,  the  Fabry-Perot 
wavelength  (828  nm)  is  shorter  than  that  of  the  gain  (~843  nm)  or  absorber  quantum  well  (~  851 
nm).  In  this  case,  changing  the  absorber  DC  bias  does  not  result  in  a  significant  change  in  the 
absorption  or  the  modulation  response,  as  shown  in  Figure  1 0. 

We  expect  that  ahgning  the  Fabry-Perot  wavelaigth  to  be  redder  than  the  absorber  quantum-well 
band  edge  will  result  in  excellent  modulation  depth  with  this  absorber  modulation  technique.  A 
reverse  bias  across  the  absorber  can  then  be  used  to  sweqp  the  absorption  edge  across  the  laser 
emission,  resulting  in  a  very  hi^  on-off  ratio.  Future  work  will  include  optimization  of  the 
bandwidth,  growth  and  fabrication  for  large-signal  modulation,  and  single-mode  operation  to 
demonstrate  low  frequency  chirp. 

1.7.  Conclusions 

In  conclusion,  we  have  constructed  a  novel  VCSEL  with  a  voltage-controUed,  quantum-well 
absorber  int^ated  into  the  mirror  stack.  If  the  device  is  designed  so  that  the  absorber  wavelength 
is  sh^tly  shorter  than  the  Fabry-Perot  wavelength  under  no  appUed  bias,  the  device  exhibits  a 
region  of  sharp  negative  differential  resistance  in  the  absorber  current  versus  voltage  characteristic. 
By  choosing  the  absorber  bias  conditions  so  that  the  operating  load  line  either  thrice  intersects  or  is 
tangential  to  the  absorber  k-Va  trace,  the  device  can  be  made  to  exhibit  either  optical  bistability  or 
self-pulsation,  respectively. 

Using  this  device,  we  produced  for  the  first  time  tunable  self-pulsations  with  an  extremely  narrow 
RF  Unewidth  of  700  kHz.  By  simply  adjusting  the  bias  conditions,  the  GHz-range  self-pulsation 
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frequency  can  be  tuned  over  700  MHz.  By  careful  device  design,  the  VCSEL  can  be  made  to  self- 
pidsate  at  low  absorber  biases,  GHz-range  frequencies,  and  potentially  at  hi^  output  power. 


Alternatively,  the  device  can  be  designed  with  the  lasing  wavelength  slighdy  shorter  than  the 
absorber  wavelength.  Using  this  device  desigi,  we  have  demonstrated  a  novel  tedmique  of 
intracavity  quantum-well  absorber  modulation  of  a  VCSEL.  We  achieved  with  this  method 
reasonable  small-signal  bandwidths  of  9  GHz  and  a  bandwidth  dependence  on  puirq)  current  of  7 
GHz/VmA.  We  analyzed  the  chirp  of  an  intracavity  modulator  in  comparison  with  that  of  direct 
modulation  and  showed  significantly  less  chirp  in  our  technique  at  hi^er  frequencies.  Future  work 
will  include  optimization  of  the  bandwidth,  growth  and  fabrication  for  large-signal  modulation,  and 
sin^e-mode  operation  to  demonstrate  low  frequency  chirp. 
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2.  In-building  Radio  Distribution  Using  a  "Saturated"  Distributed-Antenna 
Architecture 

Section  abstract 

The  radio-propagatictti  environnient  inside  buildings  is  of  sufficient  complexity  such  that  the  design 
of  an  in-building  wireless  communication  system  often  necessitates  site-specific  radio 
planning/modeling/measurements.  A  distributed-antenna  architecture  is  often  employed  for  in¬ 
building  radio  distribution.  Sophisticated  simulation  tools  have  been  developed  to  determine  the 
proper  placement  of  the  antennae.  With  the  recent  emergence  of  very  low  cost  distributed-antenna 
hardware,  it  becomes  economically  feasible  to  blanket  a  building  with  a  large  number  of  antamae 
elements  (the  "saturating  approach").  This  offers  the  possibility  of  negating  the  need  for  site- 
specific  planning,  modeling  and  survey  for  d^loyment  of  in-building  wireless  systems.  We 
investigate  the  characteristic  radio  distribution  statistics  under  this  scenario,  and  found  that  with  a 
"sufficiently  large"  number  of  antennae  elemaits  per  unit  area,  the  radio  statistics  approaches 
certain  limiting  distributions  vdiich  are  largely  independent  of  the  specific  floor-plan.  These 
limiting  distributions  can  be  accurately  predicted,  without  invokiiig  adjustable  parameters,  using 
only  a  few  key  "global"  building  parameters  vdiich  can  be  easily  determined  for  a  given  building. 
Quantitative  results  from  simulation,  as  well  as  vdiat  constitutes  a  "sufficiently  large"  number  of 
antama  elemaits,  are  described  in  detail  in  this  section. 

2.1  Introduction 

The  next  generation  of  wireless  radio  systans  promises  to  give  consumers  the  luxury  of 
tether-less  access  to  telephone  conversation,  email,  and  network  services.  It  is  clear  that  to  meet  the 
demand  of  these  new  services,  cellular  and  PCS  service  providers  must  i^jgrade  their  networks  to 
provide  complete  radio  coverage,  both  outdoor  and  indoor.  Outdoor  wireless  coverage  is  currently 
provided  by  macroceUular  base-stations,  siq>plemented  by  coaxial  or  fiber-optic  analog  transport 
links  (the  wired  RF  transport  infrastructure)  wbich  distribute  the  RF  signal  to  areas  of  poor 
coverage  such  as  subway  tunnels,  shopping  malls  and  shadowed  r^ons[l-3].  In  addition  to 
providing  a  more  uniform  radio  coverage,  this  type  of  "distributed-antenna"  network  can  be  used  to 
optimize  radio  Channel  utilization  by  dynamically  allocating  chaimel  resources  according  to  local 
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traffic  pattern  throu^out  the  course  of  the  day. 


The  concept  of  applying  distributed-anterma  architecture  to  in-building  radio  distribution  has  been 
proposed  and  discussed  in  considerable  detail[4,5].  For  practical  deployment,  it  is  necessary  to 
address  the  following  issues.  First,  the  fiber-optic  and  coaxial  cable  hnks  used  for  outdoor  coverage 
are  hi^  performance,  hi^  dynamic-range  systems  and  thus  are  expensive,  therefore  it  is  not 
economically  feasible  to  deploy  a  large  number  of  such  hnks  and  antennae  inside  a  typical  office 
building.  On  the  other  hand,  the  in-building  radio  propagation  environment  is  a  very  conqjlex  one, 
necessitating  radio-modeling/plarming  using  sophisticated  simulation  tools,  and/or  site  survey,  in 
order  to  decide  on  the  proper  locations  of  the  antennae.  This  is  exacerbated  by  the  hi^  cost  of 
instaUing  single-mode  fibers  and  heavy  coaxial  cables  in  odd  locations  inside  buildings  (as  dictated 
by  radio  propagation  requirements)  where  there  are  no  existing  conduits  or  cable  troughs.  AH  these 
issues  combined,  raises  substantial  obstacles  and  challenges  for  realizing  truly  ubiquitous  and 
seamless  wireless  coverage  of  the  next-generation  cellular  networks. 

2.2.  "Saturated"  distributed-anteima  architecture 

Instead  of  minimizing  the  number  of  antenna  deployed  in  a  building  to  reduce  cost,  we  propose 
a  contrary  paradigm  of  using  a  large  number  of  antainae  for  in-building  coverage.  This  becomes 
practically  feasible  with  the  recent  availabihty  of  very-low  cost  distributed-antenna  hardware[6]. 
We  argue  that  this  approach  will  minimize  the  total  cost  of  deploying  the  system  instead  of 
increasing  it,  as  follow: 

1 .  The  overall  radio  power  distribution  inside  the  building  becomes  more  uniform  as  the  density  of 
antenna  is  increased,  thus  reducing  the  dynamic  range  requirements  of  both  the  individual 
antenna  units,  and  the  analog  cable  links  that  interconnect  them.  (Another  way  to  state  this  is 
that  as  the  area  covered  by  each  antaina  is  reduced,  the  performance  requirements  are 
correspondingly  reduced.)  This  translates  into  a  cost  reduction  in  hardware,  particularly  since 
single-mode  fibers  and  heavy  coaxial  cables  are  no  longer  needed  in  order  to  meet  the  stringent 
RF  transport  requirements  otherwise. 

2.  While  detailed  radio  planning  is  needed  to  cover  a  building  with  a  small  number  of  antennae 
due  to  large  variations  between  building  structures  and  floor  plans,  we  found  that  the  radio 
distribution  becomes  independent  of  the  exact  floor  plan  and  the  placement  of  the  antennae. 
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when  a  "large"  number  of  antenna  is  deployed.  Under  this  circumstances,  we  found  it  suffice  to 
specify  only  a  few  ^obal  parameters  that  characterize  the  building  vdiereby  the  radio 
distribution  statistics  of  the  entire  building  can  then  be  predicted  with  a  hi^  level  of 
confidence.  We  will  substantiate  this  assertion  in  the  sections  to  follow.  This  constitutes  the 
bulk  of  the  content  of  this  section. 

3.  Under  this  scenario,  there  should  be  no  need  to  perform  detailed  radio  modeling  or  planning  in 
order  to  deploy  an  in-biulding  distributed-antenna  system.  The  ^obal  building  parameters 
yields  the  global  radio  statistics  and  die  dynamic  range  performance  required  of  the  radio 
hardware,  as  well  as  the  number  of  antennae  needed;  deployment  of  the  antennae  and  the  wired 
RF  transport  infrastructure  can  be  carried  out  in  existing  building  conduits  and  cable  trouts 
since  the  exact  locations  of  the  antennae  now  becomes  non-critical,  the  ease  of  deployment  and 
the  reduced  hardware  requirements  contribute  to  reduced  installation  cost  for  the  whole  system. 

2.3  In-building  radio  modeling  and  simulation 

Previous  work  on  in-building  radio  modeling  involved,  in  one  approadi,  establishing  a 
power-decay  law  of  the  path  loss  vdiereby  the  relation  between  the  power-decay  exponent  and  the 
gaieral  building  type  is  estabUshed  by  experimental  measurements[7,8].  The  radio  power  at  a  given 
point  in  the  building  is  determined  simply  by  measuring  its  distance  to  an  antenna.  For  a  distributed- 
antenna  system,  one  simply  sums  the  radio  powers  emitted  by  the  multiple  antamae  using  the 
power-decay  law.  An  in^iroved  version  of  this  approach  includes  wall  attenuation  by  counting  the 
number  of  walls  separating  the  point  of  observation  and  the  antennae  sources.  Fading  statistics  due 
to  multiple  antenna  sources  or  reflection  from  walls  are  not  taken  into  account,  viuch  can  be  a 
serious  deficiaicy  in  modeling  a  distributed-anterma  system.  A  second  apiproach  uses  ray-tradng 
simulation  algorithm  to  predict  the  radio  distribution  based  on  the  building  floor  plan,  which 
includes  such  detailed  information  as  locations  and  reflection  properties  of  all  walls  and 
partitions[9-ll].  These  approaches  are  deterministic  and  are  appropriate  to  studying  site-specific 
radio  propagation  effects  in  a  given  building  served  by  a  small  number  of  antennae,  of  which  the 
detailed  floor  plan  and  radio  characteristics  of  building  materials  are  knowffi  to  a  hi^  d^ree  of 
certainty.  Sophisticated  simidation  tools  have  already  been  developed  for  performing  these  planning 
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tasks. 


An  alternative  approach  is  to  deal  with  the  radio  coverage  problem  on  a  statistical  basis.  Instead  of 
attempting  to  predict  the  radio  power  at  a  specific  point  in  space  inside  the  building,  the  statistics 
over  the  entire  building  is  expressed  as  a  distribution,  from  vsfuch  the  ^obal  performance  of  the 
wireless  communication  system  such  as  dynamic  range  requirements  and  call-blocking  probability, 
etc.,  can  be  predicted.  It  is  a  useful  approach,  particular  since  intense  computation  is  not  required  as 
in  ray-tracing  approach,  but  only  if  such  radio  statistical  distributions  can  be  predicted  with 
sufficient  accuracy.  These  statistical  distributions  are  commonly  derived  based  on  ideal  models  of 
physical  processes  (such  as  addition  of  random  vectors  of  equal  magnitudes  for  Raylei]^ 
distribution,  that  with  a  dominant  component  for  Ridan,  and  others  such  as  Nakagami,  etc., [12- 
16]).  As  in  the  power-decay  law  approach,  these  ideal  physical  models  are  too  simplistic  to 
incorporate  actual  building  parameters  in  a  meaningful  way,  and  as  such  does  not  reaUy  provide  an 
alternative  to  detailed  ray-tracing  for  actual  design  purpose. 

2.4  Statistical  RF  modeling  of  in-boilclii^  distributed-antenna  system 
2.4.1  Statistical  versus  "exact**  approaches 

A  prindple  goal  for  site-specific  radio  simulations  is  to  determine  proper  placement  of 
antennae  throughout  the  building  so  that  adequate  radio  coverage  is  guaranteed  over  virtually  aU 
locations  within  the  building.  Intuitively,  such  guarantee  can  be  given  without  performing  any 
detailed  simulations  if  a  large  number  of  antennae  can  be  employed  to  "saturate"  the  budding,  radio¬ 
wise.  Performing  radio  simulation  of  such  a  scenario  using,  for  example,  WiSE  tool[10],  will  yield  a 
radio  power  distribution  throu^out  the  building  with  a  probabdity  density  function  (PDF)  vdnch 
approaches  a  5-fiinction,  as  the  number  of  antennae  per  unit  area  increases  indefinitely.  This  is  not 
observed  in  practice  since  fast-fading  due  to  multipath,  and  field  interference  between  nei^boring 
antennae,  are  not  taken  into  account  in  the  simulation.  In  fact,  we  wdl  shown  below  that  the  radio 
power  PDF  will  asymptotically  approaches  a  distribution  form  with  finite  width.  The  form  of  the 
distribution  is  largely  independent  of  the  detaded  floor  plan  but  rather,  depends  only  on  a  few 
statistical  budding  parameters  which  can  be  easdy  detennined  from  a  given  floor  plan. 
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Thus,  using  the  "saturated"  design  paradigm  for  in-building  distributed-antama  system,  the  design 
focus  should  move  from  detailed  site-specific  radio  simulation  towards  understanding  and  applying 
these  limiting  PDF's  in  the  radio  statistics.  The  rest  of  this  section  is  devoted  to  this  latter  subject. 

We  limit  our  considerations  to  the  case  \\bere  all  mobile  users  and  antennae  elements  in  the  building 
receive  and  transmit  in  the  same  radio  band,  i.e.,  the  entire  building  constitutes  a  single  "ceU"  in  the 
radio  sense.  It  is  assumed  that  RF  signals  transmitted  by  mobile  users  inside  the  building  are 
received  simultaneously  by  multiple  antenna  elements,  which  are  transmitted  via  analog  cable  or 
fiber-optic  finks  to  a  RF  summing  node  where  the  signals  from  aU  antennae  elements  inside  the 
building  are  vector-summed.  For  the  downlink,  the  RF  signals  are  distributed  from  the  central  node 
to  aU  antenna  elements  and  broadcasted  simultaneously.  As  such,  the  RF  power  distribution  so 
derived  applies  to  the  site  of  the  mobile  user  for  the  downfink  and  well  as  i^fink.  If  there  are 
inteUigait  signal  combining  techniques  at  summing  node  such  as  those  discussed  in  [17],  then  the 
analysis  should  only  applied  to  the  downlink,  and  alternative  analysis  should  be  performed  for  the 
uplink. 

Extensive  simulations,  as  described  below,  show  that: 

1 .  For  a  distributed-antenna  system  in  the  "saturated"  limit  (the  number  of  antenna  per  unit  area 
required  to  approach  this  limit  will  be  described  later),  the  resulting  radio  statistical 
distribution  is  largely  independent  of  the  exact  floor  plan. 

2.  Instead,  the  "saturated"  radio  distributions  can  be  accurately  predicted  from  these  simple 
averaged  building  parameters:  the  mean  {d)  and  variance  (a^)  of  the  distances  between  walls, 
the  average  radio  power  reflection  (R)  and  transmission  (7)  coefficients  of  the  walls. 

The  wall  separation  parameters  and  can  be  determined  quite  easily  given  the  building  floor  plan. 
The  wall  transmission  and  reflection  characteristics  can  be  found  from  well  documaited  sources  if 
the  wall  material  is  known,  and  otherwise  can  be  determined  quite  easily  throu^  simple 
measurements.  We  found  that  this  surprisingly  oversimplified  approach  can  produce  accurate 
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results  for  a  wide  range  of  building  parameters.  We  will  discuss  circumstances  under  vdiich  this 
simplified  approach  ceases  to  be  accurate. 

2.4.2  Equivalent-grid  Model 

We  demonstrate  the  validity  of  the  above  assertions  by  taking  a  typical  office  floor  plan  (for 
illustration,  that  of  Cory  Hall  at  U.  C.  Berkeley)  and  calculate  the  radio  statistical  distribution  by  (1) 
performing  an  exact  ray-tradng  calculation  on  a  floor  covered  by  a  distributed-antenna  system  with 
various  number  of  antennae,  and  (2)  perform  the  same  calculation  on  a  fictitious  floor  plan  which 
consists  of  a  rectangular  grid  of  walls  (Fig.  1 ).  The  rectangular  grid  is  uniformly  spaced  with  average 
X-  and  y-  separations  conforming  to  the  corresponding  mean  wall-sqjarations  of  the  actual  floor  plan. 


Floor  plan,  Coty  Hall  Equivalent  Grid 


Fig.l2Transformation  from  original  floor  plan  to  equivalent  grid. 

During  the  simulation,  however,  the  position  of  the  grid  line  is  displaced  from  its  r^ular  position 
according  to  a  uniform  random  variable  whose  variance  also  conforms  to  that  of  the  actual  floor  plan. 
The  actual  building  is  assumed  to  be  constructed  from  materials  with  RP  reflection  and  transmission 
coeffidents  of  R  and  T,  the  same  parameters  are  used  in  the  fictitious  rectangular  grid.  We  note  that 
even  thou^  exact  ray-tradng  are  done  both  on  the  original  floor  plan  and  the  fictitious  one,  the 
computing  time  required  to  complete  the  simulation  for  the  latter  is  about  an  order-of-magnitude  less 
than  that  of  tiie  former,  due  primarily  to  the  extreme  ease  of  determining  the  number  of  walls 
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traversed  by  each  ray,  and  locations  of  image  antennae  to  account  for  reflection  of  rays.  (We  included 
betweei  4-8  images  (reflections)  which  we  found  to  be  sufficient  for  convergence  of  the  final  result 
for  most  building  parameters.)  The  computational  resource  differential  is  particularly  significant  as 
the  number  of  antama  becomes  large. 


(M.5fn  d=7m  tf=9.4m 


Fig.l3Three  floor  plans  with  different  average  wall  separation  based  on  which  ray-tracing  and  grid 

simulations  are  done  and  compared. 


The  simulation  is  done  at  a  wavelength  of  0.3m,  wfiich  corresponds  to  a  carrier  frequency  of 
1  OUz.  The  resolution  of  the  simulation  is  one  third  of  the  wavelength,  which  is  fine  aiou^  to  resolve 
fading  effects.  The  size  of  the  building  is  50x65m,  We  compare  three  floor  plans  with  mean  wall 
separations  of  9.4m,  7m  and  4.5m,  with  standard  deviations  of  6.8m,  5.9m  and  4.5m  respectively 
(see  Fig.  2).  We  also  consider  (power)  reflection  coefficients  (7?)  of  0.1,  0.3  and  0.5,  with 
transmission  coefficients  {T=l-R).  We  explore  the  effect  of  deploying  5  to  20  antennae  which  are 
placed  approximately  uniformly  across  the  floor  as  shown  (Fig.  1).  It  is  assumed  that  the  antennae 
are  placed  on  flie  ceiling  2m  above  the  observation  plane.  The  simulation  is  a  2-dimaisional 
propagation  simulation  based  on  the  vector-sum  of  fields  at  the  observation  point. 

P  =  1  Eoj  exp(ikr,)lri  f 

where  P  denotes  the  total  power  at  the  observation  point,  n  denotes  the  distance  the  /th 
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ray  travelling  from  the  antenna/image  to  the  observer  (accounting  for  the  2m  vertical 
separation),  a,  denotes  the  attenuation  factor  due  to  the  transmission  through  walls  for  the 
rth  ray.  The  total  RF  power  emitted  from  ALL  antennae  is  assumed  to  be  a  constant 
regardless  of  the  number  of  antennae,  and  is  uniformly  distributed  among  all  antennae. 

2.5.  Results 

The  principle  results  of  the  siniulations  are  shown  in  Fig.  3  and  4.  Figure  3(a)  shows  the 
comparison  between  the  RF  power  distribution  statistics  of  the  actual  floor  plan  and  that  of  the 
corresponding  rectangular  grid,  for  three  different  floor  plans  (Fig.  2)  with  average  wall 
separations  of  9.4m,  7m  and  4.5m,  and  for  T  ranging  from  0.9  to  0.5.  Five  antennae  placed 
according  to  Fig.  1  are  used  in  this  simulation.  The  fact  that  it  is  difficult  to  differentiate  the 
existence  of  two  curves  in  most  of  the  plots  testify  to  the  surprising  accuracy  of  the  grid  simulation 
as  compared  to  that  of  the  real  floor  plan,  for  the  wide  ranges  of  building  parameters  used.  To 
further  illustrate  the  degree  of  match  between  the  two  simulations,  the  same  residts  are  plotted  as  a 
cumulative  distribution  function  (CDF)  in  Fig.  3(b)  with  a  logarithm  vertical  scale.  One  observes  a 
superb  match  down  to  a  cumulative  statistics  of  0. 1  %  at  the  tail  end  of  the  distribution.  Four 
observations  are  worth  noting  here: 

1 .  The  "smoothness"  of  the  distribution  curve  is  somewhat  surprising,  for  both  the  simulation 
involving  the  actual  floor  plan  and  the  rectangular  grid,  but  perhaps  more  so  for  the  former,  since 
the  actual  floor  plan  is  quite  irregular  (Fig.  2),  and  the  resulting  radio  power  distribution  should 
also  be  expected  to  behave  as  such.  This  "smoothness"  is  largely  a  consequence  of  a  "large" 
number  of  antenna:  if  the  number  of  anteima  is  reduced  to  below  5,  irregular  spikes  can  be 
observed  in  the  radio  power  distribution.  In  fact,  if  a  sin^e  antenna  is  used  to  cover  the  entire 
floor,  not  only  is  the  width  of  the  distribution  widened  considerable,  the  distribution  itself  is  very 
irr^ular  as  expected  of  a  conqjhcated  internal  building  structure 
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Ilg.l5(a)  Radio  statistics  (PDF)  comparison  for  two  floor  plans  with  different  antenna 
numbers  -actual  floor  plan  vs  equivalent  grid  model,  (b)  Radio  statistics  (CDF) 
comparison  for  two  floor  plans  with  different  antenna  numbers  -actual  floor  plan  vs 
equivalent  grid  model. 


2.  For  system  design  purpose,  it  is  important  to  determine  the  dynamic  range  D  of  the  received 
RF  power  with  a  given  probability  q,  i.e.,  for  a  user  randomly  located  inside  the  building,  the 
probabiUty  being  q  for  the  maximum/minimum  ratio  in  the  received  RF  power  to  be  less  than 
D  (measured  typically  in  decibel  differaice).  If  the  hardware  comprising  the  transmit,  receive 
and  transport  infrastructure  of  wireless  system  are  designed  to  have  a  dynamic  range 
performance  of  A  then  the  probabihty  of  failure  is  b-l-q,  vbich  can  be  interpreted  loosely  as 
call-blockage  probability  or  grade-qf-service  of  the  in-building  system.  Observing  the 
distribution  plots  of  Fig.  2,  one  can  discern  a  clear  trend  towards  a  narrower  distribution  (i.e.. 
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a  smaller  dynamic  range)  for  a  building  with  larger  rooms  and  lower  reflections  from  walls. 
This  can  be  intuitively  understood  since  both  of  these  trends  tend  to  diminish  effects  of  radio 
non-uniformity  arising  from  the  internal  structure  of  the  building. 

3 .  Thus  intuitively,  one  can  expect  an  empty  building  with  no  internal  structure  at  all  to  have  the 
narrowest  RF  power  distribution,  in  fact,  in  the  limit  of  an  infinite  number  of  antennae,  a  8- 
functian  for  the  PDF.  This  is  not  the  case  in  our  simulation,  as  illustrated  by  the  simulation 
results  of  Fig.  3  which  shows  the  PDF  and  CDF  of  the  radio  power  distributions  with  the 
number  of  antenna/floor  increasing  from  5  to  20,  for  two  different  sets  of  building  parameters. 
Again,  the  match  between  the  simulation  of  the  actual  floor  plan  and  that  of  the  equivalent 
grid  is  inq)eccable,  but  the  more  interesting  observation  is  the  limiting  distribution  to  which 
the  PDF(CDF)  converges  in  the  limit  of  a  large  number  of  antennae  and  an  empty  building 
(the  case  of  20  antennae  and  9.4m  rooms  with  reflection  of  0.1).  Further  increases  in  the 
number  of  antenna  and  reduced  internal  budding  structure  does  not  lead  to  significant  changes 
in  the  distribution.  This  is  due  to  radio  interference  between  antennae,  which  produces  fading 
of  the  RF  signal  in  the  space  domain,  hi  fact,  we  show  in  Section  3  that  the  distribution  in  this 
limit  actually  approaches  a  Rayleigh  distribution.  This  is  the  narrowest  distribution  one  can 
expect  for  an  m-budding  environment;  any  addition  of  internal  structures  results  in  broadening 
the  distribution  from  the  Rayleigh  limit. 

4.  For  the  grid  simulation,  the  "smoothness"  of  the  radio  power  distribution  is  seen  only  by 
incorporating  the  random  dithering  in  the  grid  positions,  as  discussed  in  Section  4  above.  The 
irr^ularity  resulting  from  stmxdation  of  a  perfectly  rectangular  grid  without  random  dithering 
is  shown  in  Fig.  5,  and  is  believed  to  result  from  the  periodic  nature  of  the  reflectors.  (As  one 
can  see,  the  effect  is  more  obvious  for  stronger  reflection  case.)  However,  the  distribution  is 
not  sensitive  to  the  standard  deviation  of  the  random  dither,  as  long  as  it  exceeds  some  value, 
typically  in  the  l-2m  range.  Therefore,  in  the  global  characterization  of  the  building,  the 
variance  of  the  wall  separation  is  not  as  important  as  the  other  three  parameters  (mean 
distance  betw/eai  walls  d,  reflection  and  transmission  coefficients  R  and  T.) 
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Fig.  16  PDF  obtained  from  ray-tracing  of  actual  floor-plan,  and  grid  model  without  random 
dithering  (simulation  done  on  the  2nd  floor  of  Cory  Hall,  Fig.  1.)  Left:  R  =10%,  7=90%; 
right:  R  =50%,  7=50% 


cN5.8m  £#=9.5m 


Fig.  17  Two  very  non-uniform  floor  plan  with  average  wall  separations  of  5.8m  and  9.5m 

5.  We  have  explored  vdiether  or  not  the  grid  model  apphes  to  an  actual  floor-plan  vdiich  is 
highly  non-uniform,  i  .e.,  there  is  a  big  disparity  in  room  sizes  in  different  parts  of  the  building 
such  as  those  with  open  atriums  and  conferaice  halls,  as  shown  in  the  two  floor-plans  in  Fig. 
6.  Results  of  ray-tracing  simulaticms  of  these  actual  floor-plans  is  compared  to  the  equivalent 
grid  model  is  shown  in  Fig.  7  (using  5  antennae  placed  in  the  same  position  as  that  before,  see 
Fig.  1)  We  observe  that,  again,  the  matching  is  very  good  for  wall  reflection  coefficients  as 
hi^  as  0.5. 


29 


PDF 


R=0.1.  7^=0.9 


R=0.3,  r=0.7 


R=0.5,  7=0.5 


«=0.3.  T=0.7  R=0.5,  T=0.5 


ijlliajiljiiijjgiliiqiili: 

. » . . . j . 

- .K— — 4^- — j - j - }•— — 

t  "  "i  i  i  i 

Fig.  18  (a)  Radio  statistics  (PDF)  comparison  for  two  very  non-uniform  floor  plans  as 
shown  in  Fig.  6  with  various  wall  reflections  and  transmissions  -  actual  floor  plan 
simulation  vs  equivalent  grid  model,  (b)  Corresponding  CDF's  of  those  shown  in  (a) 


2.6.  Further  discussions 

From  a  system  viewpoint,  the  significant  aspect  of  the  radio  power  statistical  distribution 
hes  in  the  low-power  ti^il  end  of  the  distribution,  since  these  deep  fades  corresponds  to  locations  in 
the  building  where  the  received  radio  power  falls  outside  the  dynamic  range  winch  the  hardware 
rstn  handle.  It  inqiacts  directly  the  grade-of-service  of  the  wireless  network.  To  achieve  a  1% 
grade-of-service  (i.e.,  probabihty  of  blocked/dropped  calls  is  below  1%)  which  is  ordinarily 
considered  quite  satisfactory  for  wireless  commumcation,  the  probabihty  of  occurrence  of  such 
deq)  fades  should  be  kept  weU  below  1%.  We  have  shown  that  increasing  the  number  of  antenna 
can,  15)  to  a  certain  point,  reduce  the  width  of  the  distribution  and  thus  reduces  the  probabihty  of 
deep  fades.  These  results  are  presented  again  in  Fig.  8  in  a  shghtly  differaitly  format,  in  which  the 
probability  of  the  radio  power  falling  below  a  giv®i  factor  from  the  maximum  of  the  distribution,  is 
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plotted  against  the  number  of  antennae  deployed  (for  a  50x65m  building). 


Fig.l9Percentage  of  tail  area  of  PDF  as  a  function  of  number  of  antennae  for  different 
thresholds,  (for  2nd  floor  of  Cory  Hall,  </=7m,  R=3(fi/o,  T=7(fi/o) 


S  antennae 


25  antennae 


Fig.  20  Radio  map  of  locations  where  power  i 
second  floor  of  Cory  Hall.  (rf=7m,  R=30%, 
under  "deep-fade",  2%,  right;  25  antennae 


len) 

1 35dB  below  maximum  ("deep-fade")  for  the 
70%)  Left:  5  antennae  case,  percentage  of  area 
e,  percentage  0.4%. 


From  a  mobile  user’s  viewpoint,  the  ^atial  distribution  of  deep  fades,  in  addition  to  the 
probability  of  occurrence,  also  plays  a  role  in  determining  the  quality  of  service.  Extended  zones  of 
poor  radio  perception  tends  to  be  more  annoying  than  occasional  bursts  of  noisy  pickiq)  vdien  the 
user  passes  throu^  a  deep  fade  of  limited  spatial  extent.  To  illustrate  this  aspect  of  radio 
reception,  the  ^atial  distribution  of  deep  fades,  defined  here  as  locations  wfiere  the  received  fidd 
strength  drops  bdow  35dB  from  the  ^obal  maximum,  is  plotted  in  Fig.  9  for  a  5-antenna  and  a  25- 
antenna  system.  For  the  5-antenna  system,  deep  fades  occur  over  approximatdy  2%  of  floor  i^ace, 
iduch  is  reduced  to  0.4%  for  a  20-antenna  system.  From  these  plots  one  can  notice  that  the  deep 
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fades  are  more  scattered  for  the  25-antenna  system. 

To  further  characterize  the  spatial  distribution  of  deep  fades,  we  plot  in  Fig.  10  the  2- 
dimensional  spatial  autocorrelation  of  the  deep  fades  with  a  resolution  of  1cm.  These  plots  show 
that  employing  a  larger  number  of  antennae  (approximately  15-25  antenna  over  a  50x65m  floor 
space)  reduces  the  spatial  extent  of  most  deep  fades  to  <4cm.  These  results  illustrate  the  fact  that, 
since  deqs  fades  are  caused  by  near-complete  cancellation  of  field  vectors  at  the  observation  point 
due  to  various  antenna  sources  and  multipaths,  the  spatial  extent  over  vdiich  it  occurs  is  reduced  as 
the  number  of  field  vectors  contributing  to  the  sum  increases. 


4  antennae  25  antennae 


Fig.  21  2-D  spatial  autocorrelation  of  deep-fade  pattern  (defined  as  locations 

where  RF  power  is  30  dB  below  maximum)  for  a  30m  by  30m 

2.7.  Conclusion 

In  order  to  provide  uniform  radio  coverage  in  indoor  environment,  unplemeotation  of 
distributed-antenna  networks  has  been  proposed.  We  presented  a  grid  model  to  calculate  the  radio 
power  statistics  inside  a  building  vdiich  is  shown  to  matdi  that  of  the  exact  distribute  to  a  very  hi^ 
degree  of  accuracy.  In  this  model,  four  parameters  -  reflection  coefficient  and  transmission 
coeffidmt  of  walls,  and  average  wall  distance  and  its  variance,  provide  a  statistical  description  for 
the  vsAiole  building.  The  first  two  parameters  can  be  determined  from  the  building  types,  and  the 
last  two  parameters  can  be  determined  from  the  floor  plan.  There  is  no  adjustable  parameter  in 
our  model.  The  resultant  grid  model  can  be  simulated  very  efficiently  to  give  a  statistical 
distiibutian  for  the  building.  The  advantage  of  using  a  distributed-antenna  network  to  reduce  the 
width  of  the  PDF  of  the  distribution  and  thus  inqiroving  system  performance  can  be  observed 
readdy,  especially  in  a  hi^y  attenuated  environment.  A  further  advantage  of  enqiloying  a 
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distributed-antenna  system  is  the  reduction  of  the  physical  width  of  the  deep  fades  to  below  that  of 
a  typical  antatma  size  of  a  mobile  handset,  wtodh  implies  a  virtual  elimination  the  effects  of  deep 
fades  as  far  as  user  performance  is  concerned.  The  tool  we  have  developed  promises  to  greatly 
sinphify  the  design  of  indoor  wireless  system  on  the  radio  coverage  inside  conqjlex  building 
aivironment. 
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3.  In-building  Radio  Distribution  Using  a  "Saturated"  Distributed-Antenna 
Architecture  -  Analytic  Approach 

Section  abstract 

Section  2  described  properties  of  in-building  radio  distribiition  statistics  under  a  "saturated" 
distributed-antenna  scaiario.  The  independence  of  the  asynptotic  radio-distribution  statistics  on 
specific  floor  plan  laids  the  opportunity  to  derive  these  results  analytically,  thereby  adding  insi^ 
and  understanding  to  the  problem.  The  key  assumption  in  the  analytic  approach  is  that  when  the 
number  of  antama  elements  per  unit  area  becomes  large,  we  can  approximate  the  discrete  antennas 
by  a  uniform  and  continuous  blanket  of  antenna.  The  vahdity  of  the  assumption  is  established  by 
observing  the  fact  that  the  analytic  results  so  derived  conform  extremely  well  to  that  from  the 
computer  simulation  of  Part  I,  even  for  a  surprising  small  niunber  of  antennae  per  unit  area.  Again, 
the  excellent  match  is  obtained  with  no  adjustable  parameters  in  the  model. 

3.1.  Introduction 

Distributed-antennas  are  deployed  in  indoor  wireless  microceU  to  provide  xmiforra  radio 
coverage.  With  multiple  antomas  simultaneously  transmitting  (or  receiving)  RF  signals,  the  radio 
signal  received  (or  sent)  by  a  mobile  user  anywhere  in  the  microceU  is  the  vector  sum  of  the  radio 
signals  transmitted  (or  received)  by  each  spatiaUy  located  antenna.  The  phases  of  these  component 
antenna  signals  dqjaids  on  the  relative  distance  between  the  mobile  user  and  the  antennae,  and  on 
the  length  of  the  cables  interconnecting  the  antomae  and  the  summing  node.  As  a  result,  local 
fading  may  occur  due  to  the  interference  between  different  antama  signals.  While  the  global  radio 
coverage  tends  to  become  more  uniform  with  increasing  number  of  antennae,  local  power 
variations  due  to  antama  interference  wiU  become  more  significant.  We  have  shown  througji 
simulation,  in  Section  2,  that  beyond  a  certain  antenna  "density",  the  ^obal  radio  statistical 
distribution  converges  to  a  distribution  whidi  is  largely  independently  of  the  detailed  floor  plan  of 
the  building  but  can  be  predicted  with  a  hi^  degree  of  accuracy  using  only  a  few  global  building 
parameters.  The  purpose  of  Part  II  is  to  provide  some  analytical  basis  for  this  property,  and  to 
derive  throu^  analysis  the  converging  radio  distributions  obtained  by  simulations  in  Part  I.  It  was 
found  that  excellent  agreement  can  be  obtained  between  the  simulation  results  and  the  analytical 
results,  over  a  large  range  of  building  parameters,  and  again  with  no  adjustable  parameters  in  the 


35 


model.  These  approadies  and  results  could  lead  to  a  better  basic  understanding  of  the  problem  of 
radio  propagation  in  a  complex,  in-building  environment. 


To  grasp  the  essence  of  the  problem  without  unnecessary  comphcations,  we  make  the 
following  key  assunptions: 

1.  We  approximate  the  discrete  antennae  in  a  distributed  antenna  system  by  a  continuous 
blanket  of  antennas  of  uniform  density.  The  total  number  of  antennae  in  the  system  is  the 
density  times  the  area  covered  by  the  \\4iole  distributed  antenna  system,  nominally,  one  floor 
of  a  building  or  parts  thereof  (the  "cell").  This  assumption  is  clearly  vahd  vdien  the  total 
number  of  antennas  inside  the  cell  is  very  big.  To  see  the  limit  of  vahdity  of  the  assumption, 
the  result  so  derived  wtU  be  diecked  against  actual  simulation  of  the  floor  plan  with  discrete 
autffnna  demoits.  As  reported  below,  a  suiprisin^y  good  match  was  observed  for  as  httle  as 
5  antamae  covering  a  floor  space  of  50x50m^. 

2.  We  assume  the  cell  is  infinitely  extended  and  the  cell  boundary  condition  is  thus  of  no 
concern.  In  practice,  when  the  radio  power  decays  rapidly  from  its  point  of  origin  due  to 
partitions  and  scatterers  as  typical  in  an  indoor  wireless  environment,  only  those  antennae 
close  to  a  mobile  user  contributes  significantly  to  the  radio  reception.  If  this  characteristic 
distance  is  small  compare  to  the  dimension  of  the  cell,  then  from  each  user's  point  of  view  the 
cell  is  infinitely  extended. 

Other  approximations  are  more  mathematical  in  nature  and  wiU  be  mentioned  as  they  arise  in  the 
analysis. 

It  should  be  reminded  that  the  results  derived  below  apphes  to  the  scenario  where  the 
iqilink  RF  signals  received  at  the  antenna  elements  are  simply  vector-summed  at  a  central  node  in 
the  building,  and  the  downlink  RF  signals  are  similarly  broadcasted  from  the  central  node.  More 
intplligmt  combining  techniques  require  additional  consideration  for  the  uplink. 

The  section  is  organized  into  the  following  sections.  In  Section  2.2  we  introduce  some 
mathematical  background  and  consider  the  case  where  there  are  no  partitions  inside  the  building 
(an  en^ty  building).  In  Section  2.3,  we  derive  a  theoretical  model  for  the  received  signal  amphtude 
distribution  in  a  budding  partitioned  by  a  rectangular  grid  of  walls  (which  is  used  in  the 
simulations  of  Section  2.2).  We  then  show  in  Sectiwi  2.4  the  exceUent  match  between  the  result 
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using  the  theoretical  model  and  that  from  the  computer  simulation,  again,  with  no  adjustable 
parameters  in  the  models.  We  conclude  in  Section  2.5. 

3.2.  Empty  cell  model 

3.2.1  Rayleigh’s  random  flights[l] 

A  basic  mathematical  result  to  be  used  in  subsequent  analysis  concerns  the  statistical 
distribution  of  the  sum  of  n  random  vectors.  Consider  a  two-dimensional  random  vector 
x  =  (Xj,X2)  vdiere  both  Xj  and  Xj  are  random  variables.  The  corresponding  characteristic 
function  is  defined  as 

=  (1) 

vdiere  E  stands  for  expectation.  Now  consider  the  sum  of  n  independent  vectors  of  the  same 
magnitude  x.  For  each  one  of  them,  in  polar  coordinates,  i.  e.  with  the  definitions 

X,  =  xcosO.Xj  =  xsinO,^,  =  4cos0,^2  “  ^sin0,  (2) 

we  rewrite  the  characteristic  function  as 

=  (3) 

vdiere  O  is  the  random  variable  denoting  the  direction  of  the  vector.  Since  each  random  vector  is 
issued  in  a  random  direction,  O  should  be  uniformly  distributed  over  0  to  271.  The  characteristic 
fimction  of  a  single  random  vector  is  then 

(P(^)  =  — =  y„(x^)  (4) 

iTT  I 

where  Jq  is  the  Bessel  function  of  0*  order.  Since  each  random  vector  is  independent,  the 
characteristic  function  of  the  sum  of  n  random  variables  is  the  product  of  the  individual 
characteristic  functions  of  the  random  variables,  i.  e.  Jq”{x^)  . 

To  obtain  the  probability  daisity  function  (PDF)  /(Xj,X2)of  the  random  vector 
X  =  (x, ,  Xj ) ,  one  inverse  Fourier  transform  the  characteristic  function  : 

f{x„x,)  =  .  (5) 

Since  (p(^] , ^2  )  depends  only  on  ^  +^2^  >  again  use  the  polar  coordinates  with 

Xj  =  X  COS0,  X2  =  X  sin  0,^i  =  ^  cos(}),^2  =  ^  sin  (j)  (6) 

and  express  the  PDF  as 
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/(x.,x,)  =  /(X)  =  .  (7) 

Noting  that  the  / (Xj ,  Xj  )  is  a  function  of  the  magnitude  of  the  random  vector  x  =  |x|  only,  the 
PDF  of  X  =  |x|  is  thus 

p{x)  =  27af{x)  (8) 


3.2.2  Probability  density  function  of  RF  amplitude  distribution  inside  cell 

Here  we  consider  the  case  wiiere  are  no  partitions  inside  the  cell,  and  each  antenna  is  omni¬ 
directional.  Again,  we  assume  a  2-dimensional  propagation  model  vdiere  reflections  from  the  floor 

and  ceding  are  needed.  The  anqjlitude  of  the  RF  signal  is  assumed  to  decay  as  1  /  r“  vdiere  r  is 
the  distance  between  the  source  anterma  and  the  observation  point  and  a  is  the  decay  constant  that 
is  used  in  practice  to  characterize  the  power  decay  in  an  indoor  wireless  environment[2-4].  We 
require  a  >  1 .  It  is  further  assumed  that  the  cell  is  big  so  that  the  reflection  from  the  cell  boundary 
at  any  observation  point  inside  the  cell  is  negligible.  The  received  RF  signal  at  any  particular  point 
inside  the  cell  is  thus  the  sum  of  signals  from  all  the  antennas  in  the  cell. 

The  signal  at  a  particular  observation  point  may  be  expressed  as 


»  h 


=  e 


(9) 


\\liere  the  summation  runs  over  all  the  antennas  in  the  cell.  The  constant  ^  is  a  measure  of  signal 
arr^litude  emitted  from  each  individual  antenna  (assuming  all  antennae  are  idoitical),  co  is  the  RF 
carrier  frequency.  The  phase  changes  by  2n  vvhen  the  distance  changes  over  a  wavelength  of 
the  carrier,  which  is  much  smaller  than  the  dimension  of  the  cell.  Hence  for  different  observaticjn 
points,  can  be  cxnsidered  as  a  random  variable  uniformly  distributed  over  27r.  The  received 
signal  amphtude  at  the  observation  point  is  then 


(10) 


To  find  the  PDF  of  we  first  find  the  characteristic  function  of  the  random  vector  E  and  obtain 
the  PDF  by  Fourier  transform. 


We  first  cxmsider  the  cxnitribution  from  the  antennas  in  the  differential  annulus  at  a 
distance  r  from  the  observation  point  as  shown  in  Fig.  1 .  The  total  number  of  antennas  in  the  band 
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is  calculated  as  N  =  [ilnrdr  vdiere  N  is  the  number  and  ^  is  the  antenna  density.  The  resulting 
summation  of  the  field  contribution  from  these  antennae  is 


^  /i  ^ 
6 


(11) 


Since  is  a  random  variable  uniformly  distributed  over  27t,  the  summation  is  similar  to  the 
vector  sum  in  the  Raylei^’s  Random  Fh^s  as  discussed  in  Section  1.1.  The  correspondmg 
characteristic  fimction  is  Jq^  ((Alr^  )^) .  If  we  assume  is  large  so  that  »  1,  we  can  always 
make  N  an  even  number  by  adjusting  dr  so  that 


Since  the  total  field  contribution  from  all  the  antennae  in  the  cell  is  the  sum  of  contribution  from  all 
armuh,  the  characteristic  fimction  of  the  total  field  is  fiius  he  product  of  the  individual 
characteristic  fimctions: 

2  2 

=  =  ^  °  (13) 

whereinthelast  equahty  werqjlacedtherby  y=  A^fr’^ . 

It  is  shown  in  the  Appendix  that  the  integration 

(14) 

converges  to  a  n^ative  constant  -  Cq  where  Q,  >  0 .  The  characteristic  function  is  thus  rewritten 
as: 
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=  .  (15) 

The  PDF  of  received  sisnal  E  is  obtained  by  inverse  Fourier  transformation  of  the  characteristic 
function  as  described  in  Section  1 . 1 .  Using  Eq.  (?),  the  PDF  for  the  total  received  field  amphtude 
E=\E\  is  thus 

2  2 

p{E)  =  (16) 


3.2.3  PDF  for  two  special  cases 

For  power  decay  characteristic  ot=l  (corresponding  to  the  case  ^^ben  there  is  no 

propagation  loss  within  the  cell),  the  dharacteiistic  function  is  a  Gaussian  e  ,  and  the 

PDF  of  the  received  RF  amphtude  PDF  becomes 

p{E)  =  ^i^  (17) 

vvbich  is  the  Raylei^  distribution  with  the  characteristic  parameter  b  =  ‘\ti\xCqA^  .  If  a  =  2 
(power  decays  as  ),  the  characteristic  function  becomes  e  and  the  resulting  PDF  is 

p{E)=  (18) 

with  characteristic  parameter  c  =  %ixCqA  .  For  other  decay  characteristic  a,  closed  form  of  p(x) 
is  difficult  to  obtain.  These  two  limiting  cases  for  the  PDFs  are  illustrated  in  Fig.  2.  Ih  general, 
the  Rayleigh  distribution  has  a  small  variance.  This  follows  from  the  fact  that  an  inverse-square 
law  in  power  propagation  is  indicative  of  absence  of  absorbers  or  scatterers  in  the  cell  and  thus  the 
overall  power  distribution  tmds  to  be  more  uniform.  In  fact,  in  our  2-D  distiibuted-antama 
approximation,  one  can  expect  the  Raylei^  distribution  as  the  limiting  case  with  the  smallest 
possible  variance. 
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Power  (dB) 


Fig.  23  Two  limiting  cases  of  the  RF  statistical  distribution  in  the  "empty-cell"  model. 


3.2.4  Central  limit  interpretation 

The  caitral  limit  theorem  also  appUes  here.  We  consider  the  case  vsheii  the  antenna  daisity 
p  approaches  infinity  A  shght  modification  is  however  needed,  since  when  the  antenna  density 
approach  infinity^  there  exist  in  an  infinitesimally  small  area  around  an  observation  point  an 
infinite  number  of  antennas,  which  results  in  infinite  received  power.  To  resolve  this  mathematical 
difficulty,  we  apply  a  practical  consideration  to  our  model.  Since  there  exists  a  finite  spacing 
between  users  on  the  floor  and  antennae  located  on  the  ceding,  we  introduce  a  vertical  distance  h 

between  user  observation  points  and  antennae,  resulting  in  the  actual  distance  being 

v^here  r  is  the  horizontal  distance  betweei  the  user  and  the  antenna.  We  diall  also  require  the  total 

RF  radiated  power  from  all  antamae  in  the  cell  to  be  finite  and  constant.  Thus  the  normalized  RF 

amplitude  emitted  from  each  antenna  is  A  =  \l  ■yfju  . 

With  the  above  modifications,  the  characteristic  function  becomes 


(p(^)  =  e 


2 

2itM J"r<frlog| +h^ )“)|  _  (2jin4“/a)Jo 


v^ere,  in  the  last  equahty,  y  =  ^l{r^+h^  .  When  p  approaches  infinity,  we  can  approximate 

resulting  in  the  characteristic  function; 


<p{,4)  =  e 


M-xo  Jo  4/j  _  g  4(a-l)h‘‘“ 


(21) 


\^iiic3i  is  a  Gaussian  as  e?q)ected  from  the  central  limit  theorem.  The  resulting  PDF  of  received 
signal  arr^jHtude,  r^ardless  of  the  decay  characteristic  a,  is  Ra3dei^  with  parameter 


*  4(a-l)/j^“"^  ■ 

3.3.  CeD  with  Grid  WaDs 

Actual  building  floor-plans  are  con5)lex  and  irr^ular  so  that  th^  do  not  lend  themselves 
to  analytical  derivation  of  radio  propagation  characteristics.  However,  the  results  of  Part  I  show 
that  in  the  saturated  distributed  antenna  sceiaiio,  the  global  radio  statistics  do  not  depaid  on  the 
exact  floor  plan,  m  fact,  a  overly  sin^Hstic  floor-plan  consisting  of  a  rectangular  grid  yields  results 
which  are  amazingly  dose  to  that  of  the  actual  floor-plan.  This  empirical  property  brings  forth  the 
possibUity  of  analytical  derivation  of  the  asynptotic  statistical  distributions  confuted  in  Part  I  by 
analyzing  the  rectangular  grid  instead  of  the  conplex  actual  floor-plan. 

The  grid  model  is  shown  in  Fig.  3.  For  illustrative  purpose,  the  vertical  and  horizontal 
walls  in  the  figure  are  equally  spaced  with  distances  Xq  and  respectively,  hr  the  analysis, 
however,  we  don’t  require  the  walls  to  be  exactly  equally  spaced  and  Xq  and  rqjresent  the 
average  distances  between  walls. 
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Fig.  24  Model  with  grid  walls.  Paths  1  and  2  show  that  the  path  loss  is  direction-dependent 

In  the  presence  of  the  walls  in  the  cell,  the  signal  propagation  loss  from  the  source  antenna 
to  the  observation  point  wiU  depend  on  the  number  of  walls  encountered  in  the  propagation  path. 
Moreover,  the  signal  at  each  observation  point  is  con^iosed  of  not  only  the  signals  from  the  source 
antennas  but  also  the  reflections  from  differoit  walls  wbidi,  equivalaitly,  can  represented  by 
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images  antennae  generated  by  the  original  antennas  upon  reflections  from  dififeroit  walls.  A  further 
complication  is  that  for  each  observation  point,  only  a  subset  of  aU  image  antennae  is  effective. 


We  shall  perform  the  analysis  in  three  steps.  First,  we  ignore  the  unage  antennae  by 
assuming  the  reflection  is  very  small.  The  effect  of  image  antennae  is  then  included,  using  a 
simplistic  assumption  that  the  image  antennae  are  uniformly  distributed  throu^out  the  cell.  This  is 
not  rigorously  correct,  but,  curiously,  gives  a  better  fit  to  the  actual  simulation  results  than  the 
exact  approach,  which  is  to  take  the  rigorous  spatial  distribution  of  the  image  antennae  into 
account. 


3.3.1  Analysis  with  No  Images 

In  the  presence  of  grid  walls,  the  path  loss  is  directional  as  shown  in  Fig.  3.  For  Path  1  as 
drawn,  there  are  5  walls  encountered.  If  each  wall  has  an  amplitude  transmission  factor  t  =  -n/T  , 
then  the  resulting  field  amplitude  at  the  observation  point  is  At^  I  r  .\n  general,  if  the  observation 
point  is  at  (0,0)  and  the  source  antainae  are  at  (r,©)  in  polar  coordinates,  the  resulting  field  at 
the  observation  point  due  to  the  source  antennae  is 

log/( - ■+ - ) 

E  =  e  .  (22) 

It  is  assumed  here  that  the  transmission  and  reflection  coefficients  of  the  walls  do  not  depend  on  the 
incident  angle  Since  the  characteristic  function  of  the  sum  of  the  signals  from  the  antennas  in  the 

differential  area  rdQdr  is  ^  received  field  amplitude  characteristic  function  can 

be  calculated  as 


rdr\  d01ogJo(  — 

Jo  I  j 

(p(^)^e  =e  ' 


WCM(0-flO) 


'r^+h^ 


-i)'. 


(23) 


vdiere  constant 


can  be  interpreted  as  the  average  differential 


transmission  per  unit  distance  and  Gq  =  tan  *(Xq  / ^q)  and  \^here  we  stiU  assume  the  antennae  are 
at  a  distance  h  above  the  ground.  The  PDF  of  received  field  anplitude  is  then 


p(,E)  =  E[^^„(Ei)e 


00  If/ 2  KTCOdO-^Q) 

4fi\  rdr\  rf^log7o( —  ^  ^ 

Jo  Jo 


(24) 


43 


33.2  Analysis  with  a  uniform  distribution  of  image  antennae 

In  this  scaiario,  the  resulting  field  amplitude  at  observation  point  (0,0)  due  to  image 
antamas  at  (r,  0)  is 


E  = 


rjoos^l 

Iog/-(J - '-+-i - 

,  *0  >0 


(25) 


v^^ere  the  additional  factor  (conqiare  to  Equation  (21))  is  due  to  the  fact  that  the  image  antennae 
emits  an  equivalent  field  strength  equal  to  the  reflection  coefficient  of  the  wall.  If  we  assume  that 
the  walls  are  dithered  from  their  r^ular  grid  positions  randomly  (as  eiqilained  in  Part  I)  so  that  the 
image  antennas  are  phase-independent  of  the  original  antennas,  the  characteristic  function  is  then 
the  product  of  that  of  the  original  antennas  and  that  of  the  image  antennas,  i.e. 


log 


=  e 


,tr  co<^-tfo) 

JoC  ,  .  r-Dl 


4? 


+h^ 


/  ;  2-— ^>1 


dB 


(26) 


The  PDF  of  the  received  signal  ampHtude  is  then 


^  #a’/2 


M- 


-5) 


Ji_,2  «rco<«-«o) 


dB 


dl 


(27) 


Image  antennae 
Strenatb  ^  P  /  T 


Flg.25  Illustration  of  image  antenna  distribution  where  the  original 
antenna  is  uniformly  distributed,  and  only  one  reflection  from  the 
grid  walls  is  taken  into  account. 


33.3  Rigorous  treatment  of  image  antennae 

A  more  rigorous  approach  can  be  taken  if  we  limit  the  analysis  to  include  only  the  image 
antennae  resulting  from  a  sin^e  wall  reflection.  Let  the  observation  point  be  (0,0)  and  Consider  the 
imagfi  antprma  of  the  an  actual  anterma  placed  at  (r^y)  as  shown  in  Fig  4.  The  effective  image 
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antenna  density  due  to  horizontal  walls  is  the  overlay  of  half-planes  of  image  antennae  due  to 
successive  walls,  and  thus  the  density  of  image  antennae  increases  approximately  as 


yo 


(28) 


vdiere  yo  is  distance  between  horizontal  walls  and  p  is  the  real  antenna  daisity.  In  general  the 
effective  first-order  (single  reflection)  image  daisity  at  (x,  y)  is 


M— +— )• 


yo 


and  the  effective  first-order  image  antama  daisity  in  the  differential  area  n/dtfr  is  then 

r|cos^|  r|sin^|^ 


firdrd6(- 

The  characteristic  function  (25)  thus  becomes 


Jo 


(p{^)  =  e 


Jo(  , '  ,  d.M  -- 


VK+/I 


oos(e-eori'>8-'o( 


(29) 


(30) 


dd 


(31) 


where  .  The  PDF  of  received  signal  anqjhtude  is  then 


p{E)  =  E^^^,{E^)e 


4,Hj^rd.j‘ 


•ft  12  I  e 

log'./o(- 

'0  I 


4)idW+-l-f*  r^drf  cos(6-S(,)log 

Vr^+A^  I  '<)•'» 


lTax(S-fio) 
■^oC  /  I,  , 


dO) 
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(32) 


3.4.  Comparison  with  Simulation 

The  results  obtained  in  Section  3  are  checked  against  con:5)uter  simulations  of  the 
rectangular  grid  model.  Since,  in  Part  I,  we  have  estabhshed  the  statistical  agreement  between  the 
grid  model  and  the  actual  floor-plan,  agreement  of  the  analytic  results  denved  above  unphes  an 
agreemait  with  the  results  from  the  simulation  of  the  actual  floor-plan.  First,  the  PDF  expression 
(24)  derived  fi'om  the  image-less  model  is  examined  for  hi^  transmission  (90%)  as  shown  in  Fig. 
5(a)  where  the  average  wall  distance  in  both  horizontal  and  vertical  direction  is  7  meters.  An 
excellent  match  is  clearly  observed.  The  discrqjancy  betweoi  the  theorrtical  expression  (24)  and 
real  simulation  becomes  significant  when  the  wall  transmission  is  low  (or  whai  the  reflection  is 
high)  as  shown  in  Fig.  5(b)  and  5(c)  where  the  transmissions  are  70%  and  50%  respectively.  The 


45 


reason  for  this  discrepancy  is  clearly  due  to  the  n^ect  of  image  antennae  ^^Alich  cannot  be  ignored 
at  hi^  reflections. 


RFpowBr<dB) 


RFpoweitdB) 


RF  power(dB) 


Fig.  26  Radio  statistics  comparison  for  different  R  (wall  reflection  coefficient)  and  T  (wall 
transmission  coefficient)  where  the  dotted  curve  shows  the  ray-tracing  simulation  result  as 
shown  before  in  Part  I,  and  the  solid  curve  is  result  using  Equation  (24)  without  image 
consideration.  The  average  wall  distance  is  7  meters  and  total  number  of  antennas  is  5. 


On  the  other  hand,  the  PDF  of  Eq.  (26)  which  takes  into  account  of  a  single  overlaying 
layer  of  uniformly  distributed  image  antennas  shows  excellent  match  vwth  simulation  for  almost  all 
the  wall  transmissions  from  as  low  as  50%  iqj  to  90%,  and  for  different  average  wall  distances  as 
shown  in  Fig  6.  To  illustrate  the  degree  to  which  the  agreement  holds,  the  commutative  density 
function  (CDF)  of  the  corresponding  PDFs  in  Fig.  6  are  plotted  in  Fig.  7  down  to  a  probability  of 
10"^.  It  is  interesting  to  note  that,  this  excellent  agreemeit  is  obtained  for  as  httle  as  5  antenna 
covering  an  area  of  50x50m^.  As  expected,  the  agreement  only  inproves  with  a  hi^er  density  of 
antama,  as  illustrated  in  the  PDF,  and  the  corresponding  CDF  plots  in  Figures  8  and  9 
respectively. 


It  is  a  curious  fact  that  the  "rigorous"  results  of  Eq.  (32),  which  properly  accounts  for  reflections  of 
the  antennae  from  successive  walls,  do  NOT  agree  as  well  as  that  of  a  uniform  layer  of  rniage 
antennae  as  represented  by  Eq.  (26).  The  reason  of  this  discrepancy  is  not  well  understood  at  this 
point.  However,  taking  these  results  at  face  value,  the  analytical  model  of  a  uniform  layer  of  unage 
antennae  yields  such  excellent  matching  to  the  simulation  of  the  rectangular  grid  floor-plan  (and 
thus  the  actual  floor-plan,  per  results  of  Part  I)  that  it  can  be  utilized  to  predict  radio  statistics  of  a 
building  covered  with  a  "saturated"  distributed  antenna  network,  to  a  hi^  degree  of  certainty. 
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Fig.  27  (a)  Radio  statistics  (PDF)  comparison  for  different  d  (average  wall 

distance),  R  (wall  reflection  coefficient),  and  T  (wall  transmission  coefficient)  where 
the  dotted  curve  is  the  result  from  ray-tracing  simulation  as  shown  before  in  Part  I, 
and  the  solid  curve  is  the  result  using  Equation  (26).  The  total  number  of  antennas  is 
5.  (b)  CDF's  corresponding  to  the  PDFs  in  (a). 
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Fig.28(a)  Radio  statistics  (PDF)  comparison  for  different  antenna  numbers  where  d, 

R,T  'v&  similarly  defined  as  in  Fig.  6  and  where  dotted  curve  is  simulation  result  and 
solid  curve  is  the  theoretical  result  using  Equation  (26).  (b)  Corresponding  CDF's. 

3.5.  Conclusion 

A  theor^cal  model  for  received  RF  statistical  distribution  in  indoor  wireless  environment 
with  a  "saturated"  distributed-anterma  network  is  presented  in  this  section.  Starting  from  the 
Ra34d^’s  theory,  we  first  derived  the  signal  distribution  in  an  empty  cell  imder  the  assumpticm 
that  the  antennas  are  “continuously”  and  uniformly  distributed  in  the  cell  (a  truly  distributed 
antenna).  We  deduced  closed-form  eTqpressions  for  the  distribution  for  differmt  decay 
characteristics  inside  the  cell.  The  distribution  derived  conforms  to  the  central  limit  theorem  when 
the  antenna  density  approaches  infinity  We  then  extended  the  analysis  to  consider  the  case  wiiere 
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there  are  rectangularly-grided  walls  inside  tiie  cell.  We  derived  models  with  and  without 
considering  antarma  images.  The  model  derived  by  assuming  a  uniform  image  antenna  density 
conforms  very  well  to  the  conqiuter  simulation.  The  approach  taken  and  models  derived  in  the 
section  can  be  used  to  explain  some  observed  fact  in  Section  2  and  help  deepen  the  understanding 
on  the  complex  indoor  radio  coverage  problem. 

3.6  ./^pendix  :  Proof  of  the  Convergence  of  JJ 

We  prove  the  above  int^al  converges  to  a  n^ative  constant.  Since  JqCt)  never  bigger 
than  1,  it  is  obvious  that  the  above  expression  evaluates  to  a  native  value.  In  the  following,  we 
first  prove  the  integration  converges  at  y  =  0 .  We  then  show  the  int^ation  is  also  fimte  at  each 

singularity  of  log]  JqCt)!  •  Finally  we  prove  the  convergaice  of  the  integration  at  infinity. 

Consider  the  behavior  of  the  int^and  around  0.  When  y  is  very  small,  we  notice  that 

log|./oO’)|~-J’'  (Al) 

and  so  the  int^and  behaves  like 

(A2) 

If  we  assume  a  >  1  (inverse  square  law  in  RF  propagation),  then  convergence  is  guaranteed 
around  0,  since  1  -  2  /  a  is  never  less  than  - 1 . 

Next  we  prove  the  int^ation  is  finite  around  the  singularities  of  log|jo(T)|>  i-  ®  >  th® 
zeroes  of  Joiy) .  It  is  known  that  JqCt)  c^n  be  expressed  in  terms  of  its  zeroes,  i.  e., 

=  (A3) 

i=l 

\^here  is  the  ith  0  of  JQ^y) .  To  verify  the  int^al  is  finite  around  ,  we  need  only  to  verify 
Jy'^5  log]  Jo(t)|^  is  finite  v^ai  6  is  very  small,  which  using  the  above  product  expression  of  the 
Bessel  function,  can  be  written  as 

ly^)dy  =  \^^^^^o^  +  ylyi\dy-^f‘_g^o^-yly\dy  (A4) 

where  only  the  last  term  may  blow  ip.  We  then  express  the  last  term  as 
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^'‘^\o%\^-yly,\dy  ~ 2y,fJ''\ogxdx  =  (S /y,)(log(S /y,)-l)  (A5) 

which  is  finite  and  thus  ends  this  part  of  the  proof. 

Finally  we  consider  the  convergence  of  the  integration  at  infinity.  The  asymptotic 
expansion  of  the  Bessel  function  J^iy)  for  large  y  is 

Jo(T)~logV^-|logT+logcos(T-|) .  (A6) 

Consider  the  int^ation  J”  T  where  7  is  a  large  constant.  Using  the  above 

asymptotic  expansion,  we  approximate  it  by 

JJ'  log  yfYhrdy  +  ^  logjrfy  +  y  log  cos(y  -^)dy 

(A7) 

where  obviously  the  first  term  converges.  We  rewrite  the  second  term  as 

lj^^-(i+2/a)  =  _^^y-^>-iogy\"  -  JJ (A8) 

71 

which  is  also  bounded.  The  last  term  involves  the  function  logcos(>’-— )  with  periodic 

singularities.  Sunilarly  as  before,  we  can  prove  the  mt^ahon  log  cos(jy - )ay  converges 

J  Py  ~5  4 

7C 

at  any  singularity  j, ,  since  the  integrand  logcos(j-^)  at  any  singularity  behaves  hke 
log|>’->’y|  whose  integration  around  y^  is  bounded.  Now  we  can  set  an  vpper  bound  for  this 
integration  over  any  period,  i.e.,  T  log cos(>’ -—)«/>’  <  We  then  show  that  the 

Jy  4 

last  term  is  bounded  by  the  following  inequahties: 
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-(1+2/a) 


log 


n 


cos(j/--)^fy 


1=0 


<  «.r 


(A9) 


which  completes  the  whole  proof. 
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